Ground state cooling of mechanical resonators 
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We propose an application of a single Cooper pair box (Josephson qubit) for active cooling of 
nanomechanical resonators. Latest experiments with Josephson qubits demonstrated that long 
coherence time of the order of microsecond can be achieved in special symmetry points. Here we 
show that this level of coherence is sufficient to perform an analog of the well known in quantum 
optics "laser" cooling of a nanomechanical resonator capacitively coupled to the qubit. By applying 
an AC driving to the qubit or the resonator, resonators with frequency of order 100 MHz and quality 
factors higher than 10 3 can be efficiently cooled down to their ground state, while lower frequency 
resonators can be cooled down to micro-Kelvin temperatures. We also consider an alternative setup 
where DC-voltage-induced Josephson oscillations play the role of the AC driving and show that 
cooling is possible in this case as well. 



INTRODUCTION 

Recently, fabrication of nanomechanical resonators 
with fundamental frequencies in the microwave range 
(100 MHz to 1GHz) has been achieved 0. For such res- 
onators, the quantum mechanical level spacing is a few 
micro-eV, which is comparable to the lowest achievable 
cryogenic temperatures. Freezing out the mechanical de- 
grees of freedom is favorable for ultra sensitive detection 
applications Q due to reduced effects of thermal fluctu- 
ations. Even more spectacular applications can be en- 
visioned if it is possible to cool the mechanical systems 
down to their motional ground states with high probabil- 
ity. Creation of exotic non-classical states, entanglement 
with other quantum objects, e.g. spins or atoms, co- 
herent quantum information transfer between quantum 
sub-systems are just a few possibilities. However, reach- 
ing the motional ground state using conventional passive 
cooling techniques is practically unfeasible, and there- 
fore other approaches need to be explored. Fortunately, 
the cooling problem is not unique to nanomechanics; a 
similar problem has been encountered and successfully 
solved in the field of ultra-cold atoms, where by using 
active cooling approaches it was possible to quench the 
vibrational motion of atoms and reach effective nano- 
Kelvin temperatures Q • The connection to quantum op- 
tics has been recently explored in several works. Hopkins 
et al. 4] applied the quantum feedback control ideas to 
nanomechanical resonator cooling. Wilson-Rae et al. Q 
proposed an analogue of resolved sideband laser cooling 
by coupling the resonator displacement to the level spac- 
ing of an attached semiconductor quantum dot, which 
is being irradiated by red-detuned laser. An advantage 
over the feedback-based techniques is that the sideband 
cooling does not require on-the-fly analysis of the out- 
put of a nearly ideal detector. On the other hand, direct 
implementation of the Ref. approach appears rather 
challenging, from the fabrication stand point and due to 



stringent constraints on the quantum dot relaxation rate, 
which needs to be slower than the resonator frequency for 
the vibrational sidebands to be resolved. 

In this work we study an alternative realization of the 
laser- like cooling for nanomechanical resonators, where 
the role of the two-state system is played by a supercon- 
ducting qubit (Cooper Pair Box, or CPB) capacitively 
coupled to the resonator. Interaction between the qubit 
and the resonator leads to splitting of the qubit states 
into equidistant vibrational sidebands. Latest experi- 
ments with Josephson qubits showed that long coher- 
ence time of the order of microsecond can be achieved 
in special symmetry points 0. When the qubit relax- 
ation and dephasing rates are smaller than the oscilla- 
tor frequency, one reaches the resolved sideband regime, 
favorable for cooling. By tuning microwave source fre- 
quency into the first red sideband (qubit level spacing 
minus the oscillator frequency), one can ensure that the 
microwave photon absorption processes are preferentially 
accompanied by simultaneous phonon emission from the 
resonator. The cooling cycle is completed when a photon 
is spontaneously emitted at the qubit natural frequency 
into an external bath. The emitted photon is blue-shifted 
relative to the source, as it carries away one resonator 
phonon energy. The exact reverse of this process is ex- 
ponentially suppressed if the qubit level spacing is large 
relative to temperature, and hence the heating is deter- 
mined by other - non-resonant or driving-induced - pro- 
cesses which are much slower than the dominant cooling 
mechanism. This makes cooling possible. A suggestion 
to use coupled nanoresonator-qubit system for laser-like 
cooling has been first made by Irish and Schwab Q ■ Here 
we consider two implementation of this scheme that use 
different AC driving sources: (1) microwaves directly ap- 
plied to the CPB or the resonator, in the form of AC 
flux or voltage; (2) the AC Josephson effect on an aux- 
iliary Josephson junction of the qubit. The advantage 
of the second approach is that it only requires DC bias 
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FIG. 1: The system with AC driving. 

for cooling; however, we find that it is not as effective 
as the AC biasing scheme. In the opposite limit of fast 
CPB relaxation when the vibrational sidebands are not 
resolved, we demonstrate that another type of laser cool- 
ing, "Doppler" cooling jsj, can be performed down to 
the temperature defined by the quantum dot level width. 
Although ground state cooling in this regime is impos- 
sible, this technique could be practically attractive for 
noise reduction in local probe applications (e.g., AFM, 
MFM, and MRFM), where typical resonator frequencies 
are below 1 MHz. 

AC COOLING 

System. The system under consideration is shown in 
Fig. ^ A mechanical resonator (horizontal beam) is cou- 
pled to a Cooper pair box (yellow rectangle) through ca- 
pacitance C x (x), which depends on the resonator dis- 
placement x. The charge, and hence the state, of CPB 
is separately controlled by the gate voltage V g applied to 
capacitor C g . CPB is coupled to a large superconductor 
through two Josephson junctions. The SQUID geometry 
is chosen to allow for application of an external AC flux 
to the system to provide the AC driving needed for cool- 
ing. Similar systems were considered in Refs. The 
Hamiltonian of the system without dissipation reads 

H - «-<W-g.M*>' - «,<♦,, cos9 + H, , (1) 
20e 

where Q is the charge on the island, Cs = C + C g + 
C x (x) is the total capacitance of the island, and H x = 

2 2 2 

2_ _|_ mugx ^ rpj^ -(.Q-j-gj Josephson energy of the SQUID 
controlled by an external flux $ x is given by Ej($ x ) = 
2Ej cos(7T$ x /$o)) where E® is the Josephson energy of 
each of the junctions (we consider a symmetric setup) 
and $o = h/2e is the (superconducting) flux quantum. 
We assume the total gate charge Q g = C g V g + C x (x)V x 
to be close to an odd number of electron charges, i.e., 
Q g = 2e(N + 1/2) + 2eSN, where |<JJV| < 1/2. Then we 



can use a two state approximation, i.e., |f) = \Q = 2eN) 
and ||) = \Q = 2e(N + 1)) and rewrite the Hamiltonian 
using the Pauli matrices as 

H = AE c (x)5N(x)(r z +E c (x)(l+A5N 2 (x)) 

- —<7 X +H X , (2) 

where Eq(x) = e 2 /2Cj> The second term in the first 
line of Eq. depends on x and, thus, added to the os- 
cillator Hamiltonian H x , renormalizes (slightly) the os- 
cillator parameters. This term is also responsible for 
the direct coupling of the oscillator to the dissipation 
in the circuitry (see below). We, first, drop this term 
for clarity but later reintroduce it when discussing the 
direct coupling between the oscillator and the gate volt- 
age fluctuations. Assuming that fluctuations of x are 
small relative to the resonator-CPB distance d, we obtain 
C x (x) « C x -C x x/d. Thus SN(x) « SN-N X x/d, where 
N x = C x V x /2e, and E c {x) « E c + E c (C x /C^) (x/d). 
Then the Hamiltonian simplifies to 

E 

H = AE C 6N a z - -y a x + A (a f + a) a z + huj a) a , (3) 

w here A = - AE C [N x - SN(C X /C S )} (Ax/d) and Ax = 
h j (2mw ) is the amplitude of zero point motion, x = 
Ai(a' + a). To increase the coupling one usually applies 
high gate voltage V x so that N x 3> 1. Thus, approxi- 
mately 

\n-4E c N x ^. (4) 

Both gate voltages, V g and V x , as well as the exter- 
nal flux <& x fluctuate as they are provided by dissipa- 
tive sources. This makes 6n and Ej in Eq. © fluctuate. 
Moreover, in all real systems there are 1// charge and 
flux noises. The charge 1// noise can effectively be added 
to the noise of the gate charge 6N. Experimentally, 1/ / 
noise is the most severe factor limiting coherence. Long 
coherence times have been achieved in Ref. by operat- 
ing in a special point where the 1/ / noise is less harmful. 
In the special point, the total energy splitting of the qubit 
AE = y/(8E c ) 2 SN 2 + E'j is not sensitive to the fluctu- 
ations of SN and & x in the linear order. This implies 
(SN) = and dEj/d$ x = 0. The fluctuations of SN 
are, nevertheless, still there and we obtain the following 
Hamiltonian characterizing the special point 

Ej 

H = — — a x + hujQa^a + X(a f + a) a z 

- —a z +g(o) + a)X, (5) 

where X = (C e SV e + C X SV X ). By construc- 

tion the coupling constant A also contains a fluctuating 
part. We neglect this higher order effect. The fluctu- 
ations of voltages SV g and SV X are caused by external 
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impedances Z(uj) of the circuits that supply the voltages. 
The (unsymmetrized) noise power is given by (SV 2 ) = 
(SV(t)5V(0)) u = [RcZ(uj)}hLu{coth{hLu/2k B T) + l). Thus 
the term Xa z /2 represents the coupling of the system to 
the harmonic (Gaussian) electromagnetic bath. The situ- 
ation is similar to the quantum optics one where an atom 
is coupled to electromagnetic vacuum. The main differ- 
ences are that in our case the coupling is longitudinal, i.e., 
it does not cause spin flips in the natural (charge) basis, 
and the bath temperature cannot always be neglected. 
The last term in Eq. (J5J is the direct coupling between 
the bath and the oscillator. It originates from the second 
term in Eq. and g — N x Ax/d at the symmetry point. 
In the weak coupling case that we consider, the bath is 
fully characterized by its spectral function, and hence we 
do not need to include the self-Hamiltonian of the bath 
explicitly. The non-electromagnetic bath acting on the 
resonator is later introduced through a phenomenologi- 
cal quality factor Q. 

Equilibrium transition rates. To proceed, we first need 
to determine the resonator-qubit system relaxation rates 
caused by the electromagnetic environment. It is con- 
venient to perform a tt/2 rotation in the x — z plane, 
o~x <-> &z, to the eigenbasis of the qubit at the symmetry 
point. The Hamiltonian becomes 



H = — — a z + fiLjQa)a + A (a' + a) a x 
- —o- x + g{a) + a)X . 



(6) 



We choose Ej 3> ujq, A, T. This ensures that without the 
driving the spin is all the time in the ground state ||). 
The qubit relaxation rate is 



r r 



Ah 2 



TtCtqEj 



(7) 



where (X%) = 2irh 2 a g oj[coth(huj/2k B T) + 1], a g « 
C ^2 S jj^, Rq = h/Ae 2 (here we assumed independent 
5V & and 5V X with similar external impedances KeZ(uj) = 
R r~j 50 f2). The opposite excitation rate is exponentially 
suppressed. 

The direct coupling between the bath and the oscil- 
lator gives the dissipative rates between the oscillator 
states \n): T„^ n _i » g 2 (X 2 =Uo )n/h 2 and r„^„+i ps 
g 2 (X 2 = _ UJo )(n + l)/fi 2 . In addition, the oscillator can 
relax via the virtual excitations of the qubit. The cor- 
responding processes are shown in Fig. |2J The rates for 
these processes are given by 



r „, A 2 (XgUJ 

Wl ~£^ Ah 2 n ' 
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Ah' 



(8) 
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FIG. 2: Dissipative processes due to the presence of the qubit: 
a) n — ► n— 1; b) n — > n + 1. The spectra of the oscillator and 
the qubit are superimposed. 



These rates are bigger than those due to the direct cou- 
pling by a factor (2Ec/Ej) 2 , and hence the direct cou- 
pling term in the Hamiltonian (oc g) can be discarded. 
The above transition rates translate into the resonator 
quality factor caused by the coupling to the electromag- 
netic bath, 



1 



Qc 



El 



Ota 



(10) 



The charge noise at the (relatively low) frequency 
u> = loq may be dominated by the 1// contribution. 
The symmetrized correlator of this contribution has been 
studied, e.g., in echo experiments and one can as- 
sume S x (io) = ((X 2 ) + (X 2 J)/2 = E 2 /f /\u;\, where 

Ei/f = AEcy/oty] and typically a,\/f w 10~ 7 . We are 
not aware of any study of the unsymmetrized correla- 
tors. Thus, we will introduce a phenomenological tem- 
perature Tiji for the particular frequency to — ujq via 
exp[-^ /(fc B T 1// )] = (Xl U0 )/{XZ o ). Then, for the 
corresponding quality factor we obtain 



n/f 



8E 2 
(hujo) 2 



a 1/f 



tanh • 



TlOJQ 



2k B Ti/f 



(11) 



It is reasonable to assume that the effective 1// noise 
temperature is not lower than the environment temper- 
ature. For further analysis we will assume ~ T. 



Cooling by applying flux driving 

For cooling we need to drive the system out of equi- 
librium. We propose to apply an external AC flux, 
®x(t) = &x,o + D^o cos ujd t, to the qubit. As we op- 
erate at the special point where dEj/d$ x — 0, linear 
driving is impossible. In the quadratic order, we obtain 
Ej -> Ej + (TrD/2) 2 Ej{cos2uj d t + 1). This gives 

E J + n x t 

H = - a z - — a x + hujoa'a 

+ A (a 1 + a) a x — — <j z cos 2wd t , (12) 
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m m+1 m m+1 

FIG. 3: The cooling process. FIG. 4: The heating process. 



where ft — (irD/2) 2 Ej. To describe the cooling process 
we use the Floquet picture. In other words, we count 
the number of energy quanta taken from the pumping 
source. This amounts to substituting the factors e ±2uJd * 



by raising/lowering operators e ±x (e l 



1)) 



and subtracting the energy taken from the source, 2hu;m, 
from the Hamiltonian. Then the new Hilbert space of the 
problem is extended as \a}\n}\m), where \a) is the state of 
the qubit (||) or ||)), while n is the number of quanta in 
the resonator. Due to the new term in the Hamiltonian, 
—2hujdm, an unbound staircase of states appears. Each 
step of the staircase (Floquet zone) is characterized by 
a number m and spans the Hilbert space \a) \n). The 
neighboring zones are shifted by energy 2hu>d relative to 
each other. The states \a) |n) |m) are the eigenstates of 
the Hamiltonian 



H 



-(Ej + Q) a z /2 + Tiluqo) a — 2hLu^m . (13) 



The other three terms form the perturbation 

X Q 
H' = -— a x + X (a f + a) a x - — a z cosx 



(14) 



This perturbation H' causes transitions inside of each 
zone and also down and up the staircase. At finite (non- 
infinite) temperature of the bath the down transitions 
prevail and the system propagates down the staircase. 
This corresponds to the flow of energy from the driving 
source to the bath. During this flow, the driving source 
quanta of energy can be up- or down-converted in fre- 
quency by amount ujq, that is, photons with frequency 
2ti>d ± luq are emitted into the bath. The former case 
corresponds to cooling while the latter to heating. 

We first consider the up-conversion cooling process 
shown in Fig. We choose the optimal detuning so 
that 2cjd = Ej + £1 — loq and the levels |t,77, m) and 
|J,,7i— l,m+ 1) are degenerate (in what follows we use 
the units fi = 1 for brevity). Thus we have to be careful 
when calculating the rates. First, we obtain the second 
order matrix element between these two states equal to 



A %fn where 



A 



OA 
2Ei 



(15) 



Note that the two paths shown in Fig. [3] interfere con- 
structively. Then, the cooling rate depends on the re- 
lation between A y/n and T r /2. For the weak driving 
case Ay^i < T r /2 we have the cooling rate r^_ >n _ 1 » 
(4A 2 /r r ) n (the superscript fl is to emphasize that this 
rate is due to the driving). To justify this result we con- 
sider this process as tunneling from the level 
to the broadened level n — 1, m + 1). The retarded 
Green's function of the second level is then given by 
G R (w) = 1/(lo + iT T /2). Using the Golden rule (the 
adiabatic elimination technique [lOj), we obtain 



27rA^n 
4A 2 



--ImG K (cj = 0) 

7T 



r r 



(16) 



In the opposite, strong driving case, Ay^ > r r /2, there 
are coherent oscillations between the two levels. The 
appropriate description is then to say that doublets of 
new eigenstates, are formed which are split in energy by 
2Ay / n, see Fig. The doublets are defined for n > 1 
as = G n iT) =t \ n ~ lj-l))/V2 and for n — we 
have a single state ip n =o = 1 , T ) - There are four pos- 
sible transitions between each state of the doublet and 
each state of another doublet shifted as n — > n — 1 and 
m — > 771 + 1. All four rates are equal r r /4. In total this 



gives T„_ 



r r /2. To summarize, 



(4A 2 /r r ) ? 7 
r r /2 



if av^ < r r /2 

if A^/n~ > T r /2 



(17) 



Note that in the strong driving regime n denotes the 
doublet rather than the oscillator's level (sec discussion 
below) . 

In addition to cooling, the AC driving induces compet- 
ing heating processes. One, shown in Fig. is a stan- 
dard in quantum optics off-resonance process. Indeed, 
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FIG. 5: Level structure of the strong driving for ac cooling. 
Here tf± = (\n, T) ± \n - 1, i))/y/2. 




FIG. 6: Level structure of the strong driving for ac heating. 
The two shown processes interfere constructively. 



the states ||, n,m) and |J.,n+l,m+l) are connected 
in the second order by the matrix element A y/n + 1. 
These states are off-resonance (we restrict ourselves to 
the values of n such that level splitting 2ujq ^> A \Jn + 1) . 
Thus, in the case Ay/n < T r /2 we obtain r„^„ +1 w 
(A/2ojo) 2 r r (n + 1). In the strong driving case A y/n > 
r r /2 we use again the basis of doublets, Fig. El and arrive 
at r^ n+1 w (A/4w ) 2 r r (n + 1) (1 + 5 n , ). In this case 
the rates are between doublets as a whole. The factor 
(1 + S n ,o) is due to the fact that there is only one state 
in the "doublet" n = 0. To summarize, 



4wg 



r r (n + i) 



1 

i+<W) 



if Ay/E < r r /2 
if Ay/E > T r /2 ' 
(18) 



Another driving-induced heating process, not charac- 
teristic for quantum optics, is due to the fact that in solid 
state systems there is strong noise at low frequencies (1// 
noise). Thus, processes like the one shown in Fig. be- 
come relevant. This process excites the qubit with the 
rate 



4EI 



Ah 2 



(19) 



Since at u) = —ojq the noise is dominated by the 1// con- 
tribution, rate l|19f) might dominate the heating. With 
this process, the heating is also resonant and, in analogy 




m m+1 

FIG. 7: Qubit heating induced by the applied drive, 
to Eq. p7jl. still assuming that r e <g; T r , we obtain 



r c (4A 2 /r 2 )n if AVn<T r /2 
?f(l + 5„,o) if AV^>r r /2 



(20) 



Master equation. Taking into account the internal dis- 
sipation of the resonator with quality factor Q, one can || 
write down the master equation for the probability P n to 
find the resonator in the state |n). First, we analyze the 
(strong driving) limit, A > r r /2. Then, it is convenient 
to write the rate equation in terms of the probabilities 
D n to find the system in the doublet ip^ (n > 1). The 
probability Dq, then, is the for the system to be in the 
(non-degenerate) ground state ipo — |f, 0). The probabil- 
ities P n are given by 



Pn = g (A, 



D 



n+lj 



2 DnOn.O 



(21) 



The master equation reads 

D n = [D n+1 - D n (l - 6 n , )} 



8a; 2 



r r [nDn-i (1 + <5„,i) - (n + 1) D n (1 + S n ,o)} 



- r c [D B _ X (1 - 5„, )(1 + 5„,i) - A, (1 + *„,o)] 
uj (n B (uj ) + 1) 

" 5 X 
(n + i)A l+ i - (n - ^)A»(1 - <5„ ;0 ) 



w n_B(wo) 



n-- ) /)„ i.CI - f \,„)(l +<v,j.) 



- ) D B (1 + 5 n , Q ) 



(22) 



The unusual form of the second part of this master equa- 
tion is due to the structure of the matrix elements be- 
tween the different states of the doublets. 
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Multiplying Eqs. (|22|l by and summing over n we ob- 
tain 

d 



l -Y,{l~D )+ l -Y c (l + D Q ) 



8^ 1 + D °) 



"0 



up A 2 
Q 8u; 2 



w n B (w ) 



(23) 



where (n) = J2 n nD n- We rewrite the RHS of Eq. O 
in terms of the phonon occupation expectation value 



l-£> 



(24) 



We, then, obtain 

s« - 



16w 2 



r r (3 + £> ) 



(25) 



There are two main cooling regimes in (|25|l . If the 
coefficient in front of (n) is positive, i.e., if luq/Q 3> 
(A 2 /8cjq) r r , we obtain the usual cooling with the cooling 
rate slowing down with decreasing (n). As this is the rel- 
evant regime for realistic parameters we will analyze only 
this case. Interestingly, however, in the opposite case of 
very high Q, when ujo/Q <§; (A 2 /8w 2 ) T r , the cooling rate 
accelerates until (n) ~ 1. 

If (h) ^> 1 (to be checked for self-consistency) the prob- 
ability to be in the ground state is negligible, D n <C 1. 
Then 



(n) w n B (u} ) 



(r r - r c )Q 

2un 



(26) 



This regime, thus, is realized when jig(wo) > (r r — 
T c )Q/(2lu )- Clearly, for at least some cooling we need 
(r r - T c ) > (2uj )/Q. At lower temperatures, Eq. 
gives negative (n) which means that the approximation 
breaks down and (n) ~ 1 or less. Then, for an estimate 
we can use 1 — Do ~ (n). In this regime (n) « 2(n) (see 
Eq.lEHO, therefore 1 - D ~ 2(n). Thus, for T r > T c 



Q 



r r 

For n^fwo) > T c Q/luq we then obtain 

Lo n B {u) ) 



(>■>) 



(27) 



(28) 



while in the opposite case the average occupation satu- 
rates at 



E 



(29) 



From Eq. l|2.'{ll it is clear that for ns(wo) ^S> 1, the cool- 
ing is initially exponential in time, (n)t ~ n B (uJo) — 
r r Q/(2cjo) [1 — cxp(— yt)], with the decay rate deter- 
mined by the oscillator bare damping, 7 = loq/Q- Only 
when the low occupancy regime ((n) t < 1) is reached, 
the rate of exponential decay increases to 



7 = 



r r - r c 



U)0_ 

2Q ' 



(30) 



If the driving is weaker, then for some low values of n 
the following master equation holds (see Refs. 0,0) 



Pn = [A 



A 4 



Lo (n B (uo) + 1) 



[(n+l)P n+1 -nP n ] 



u n B (uj ) 

Q 



K„-i - (n + 1)P„] , (31) 



where A- = 4A 2 /r r and A+ = (A/2w ) 2 r r + 
r e (4A 2 /r 2 ). For the average occupation number (n) in 
the stationary state we obtain 



_|_ UJQ TI B (^0) 



A. 



A + + Q 



(32) 



Regimes similar to Eqs. (12811 and (|29J) are clearly identi- 
fied. However, a regime similar to Eq. 1(26|1 is not pos- 
sible within master equation 1)31(1. The expression for 
the oscillator occupancy Eq. IMl't can be naturally in- 
terpreted in terms of two independent heat baths acting 
on the resonator, one being the equilibrium environment 
at the nominal external temperature T, coupled to the 
resonator by a coupling strength 70 = luo/Q, while the 
other bath being introduced by the cooling process itself. 
The effective temperature of the latter, T* , can be de- 
fined through A-/A + = exp(ojo/T*), and the effective 
coupling strength is 7* = luq/Q* = A- — A + . Then, the 
final resonator occupancy can be re-expressed as 



70+7* 



(33) 



where n* B is the Bose distribution function at tempera- 
ture T*. Clearly, T* is the lowest possible temperature 
for a given cooling process, which is achieved for 7* ^> 70- 
The combined damping coefficient 7 — 70+7* determines 
the rate of relaxation to the new stationary state. 

Example. We consider a nanomechanical resonator 
with fundamental frequency 100 MHz (ujq = 2ir x 
100 MHz = 0.5 ficV) and quality factor Q = 10 5 . It 
is coupled to the qubit, which is characterized by the 
Josephson energy Ej 50 /ieV and Coulomb charging 
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FIG. 8: Cooling diagram for a 100 MHz resonator (detailed 
parameters in text). 



energy Eq r> 160 /xeV (corresponds to Cs w 500 aF). The 
coupling strength between the resonator and the qubit 
is determined by the mutual capacitance C x ~ 20 aF, 
and the gate voltage Vm ~ 1 V, such that n x Ri 60 
Cooper pairs (see Ref. 8]). The gap between the res- 
onator and the CPB is d w 100 nm. For these parame- 
ters, from Eq.|U the resonator-CPB coupling strength is 
A w 0.1 [ieV. Assuming that C g < C x , for the relaxation 
rate of CPB we find T r = 3 ■ 10~ 3 ^eV. 

The circuit-induced quality factor of the resonator is 
Qem ~ 4 • 10 9 , which is significantly higher than the qual- 
ity factors of typical resonators. For the 1// contribu- 
tion at Ti/f > 10a->o we obtain Qi/f > 10 7 , which is 
still higher than a typical value. Hence this modification 
of the oscillator damping can be neglected compared to 
other environmental effects. 

For cooling, we take w 1 ^eV, which corresponds 
to the modulation depth D = 0.1. We thus obtain 
A Ri 10 -3 /LteV, and hence for all n we have A y/n > 
r/2 and r^_ >n _ 1 Ri r r /2 (strong driving regime). The 
heating is indeed dominated by 1// noise, with r c = 
(n/2Ej) 2 (Ef /f /4wo) ~ 1-8 • 10~ 6 ^eV. In this regime we 
get the following results for cooling: 



n B (u>o) — 300 
1.5 • 10- 3 n B (Lo ) 
0.5 • 10~ 3 



if n B (uj ) > 300 

if 0.33 < n B (u ) < 300 

if 71b(o;o) < 0.33 



The exact numerical solution of the rate equations 
Eq. J23 is shown in Figure |H1 



Cooling by applying voltage driving 

Another way to achieve AC cooling is by applying radio 
frequency voltage bias to the gates. In Fig. apply a 
driving voltage V x = Vq cos lu^ t on the resonator and 
another driving voltage V g — — {C x /C g )V x on the CPB. 



The ac voltage V x generates resonant coupling between 
the mechanical resonator and the CPB when u ac — Ej — 
ujq, which corresponds to the first red sideband coupling 
in quantum optics. The voltage V x also generates an 
oscillating charge bias on the CPB with 5N* — C x V x /2e; 
however, it is balanced by the bias V g , which prevents 
harmful ac pumping of the CPB. 

The scheme. The Hamiltonian of the voltage driving 
scheme is 



H 



^-<J Z + hu>oa}a + AE c 5n q a x 



+A coswd t (jv + a) <7 X — —a x (34) 

where the qubit works at the optimal point of Sn g = 
to avoid the 1/ / noise. The coupling is 



A = -AE, 



C X V Q Ax 
2e d 



(35) 



Similar to the ac flux driving setup, we analyze the side- 
band cooling (ujq 3> A, T r ) in the regimes of weak driving 
and the strong driving; both can be reached in experi- 
ments. 

The analysis is essentially the same as that in the flux 
driving case. Instead of the second order matrix element 
A (see Eq. (|15fl ). we have the direct coupling A/ 2. Thus 
to obtain the cooling and the heating rates we substitute 
A — > A/2 into Eqs. ljT7J) and (JTSJ). In this scheme there 
is no driving induced contribution of the low frequency 
(1//) noise similar to l|20|) . 

Strong driving. The strong driving regime is achieved 
for Xy/n/2 ^> T r /2. As the 1/f noise does not con- 
tribute to the heating, i.e., the rates l|20|l do not appear, 
the leading heating mechanism is the off-resonance cou- 
pling (see Eq. l(T%|l ). Thus, the dynamics is described 
by the rate equation 122|) with T c = 0. We obtain the 
cooling results as follows: For ns > Q(T r /2uio) we ob- 
tain (n) Ri ns — Q(T r /2ujQ). In the intermediate regime 
Q(r r /2u/ ) > n B > Q(r r /w )(A/4cj ) 2 the result is 
(n) « uj n B /QT r . Finally, for n B < Q{T r / uo ){\/ Auj ) 2 , 
the occupation number saturates at (n) Ri (A/4u>o) 2 - 

For example, let r r = 5 • 10" 3 /^cV, A = 25 • 10~ 3 ^eV 
(with the voltage bias V ~ 50 mV ) and Q = 10 5 . We 
have (X/ Au!Q) 2 T r Ri 8 • 10 _7 ^eV. The stationary occupa- 
tion number is 



n B (oj ) - 500 
10- 3 n B (uj ) 
1.5 • 10~ 4 



if n B (u ) > 500 

if 0.15 < n B (ujo) < 500 

if n B (w ) < 0.15 



which gives (n) = 2.5 • 10~ 3 at the temperature of 20 mK. 

Weak driving. For X^/n/2 <C T r /2, the dynamics 
of the resonator is described by the rate equation l|31(l 
with A- = A 2 /r r and A+ = (X/Auj ) 2 T r . Note that 
the low frequency noise doesn't appear in this scheme 
(r c = 0), which improves the cooling efficiency. For 



8 



ns(wo) < A + Q/ojq the average occupation is given by 
(n) = A+/A-; for ns(oj ) > we obtain (n) = 

n B (w )woM- ( 3- As an example, let T r = 50 • 10 _3 ^ieV, 
A = 5 • 10~ 3 /ieV (with the bias voltage of Vq ~ 10 mV 
) and Q = 10 5 . We have A_ = 0.5 • lO^VeV and 
A + = 3 • 10~ 7 /ieV. Then, the cooling results are as fol- 
lows: 

(«\ = S 10 ~ 2n *(«*>)> n B (u ) > 0.06 (3g) 



6 • 10" 



n B (u ) < 0.06 



where at the temperature of 20 mK with ub(^o) — 5, we 
have (n) = 0.05. Thus better cooling can be achieved in 
the strong driving regime than that in the weak driving 
regime. 

Discussion. The relaxation rate T r can be adjusted 
by varying the external circuit of the CPB. For example, 
by varying the gate capacitance C g , the relaxation rate 
changes as T r <x Cg. The coupling constant A can be ad- 
justed by varying the bias Vq. In this scheme, we choose 
Vo to be in the range of 10 . . . 100 mV and T r in the range 
of (1 . . . 100) • 10 _3 /j,eV, which includes both the weak 
driving regime (A v / n/2 <C T r /2) and the strong driving 
regime (Ay / n/2 ^> T r /2). These parameter regimes have 
been realized in charge qubit experiments. 

One practical issue of this scheme concerns the ac- 
curacy of the gate compensation. Both gate voltages 
generate an extra part of the CPB's charging energy 
5H = AE c Sn g a x , with 5n g = (C X V X + C g V g ) /2e. By 
controlling the voltage with an accuracy of microvolts, 
which can be achieved with standard technology, the os- 
cillating bias on the CPB can be neglected. 

Compared with the ac flux driving cooling scheme in 
this paper, the resonator can now be cooled to a lower 
temperature because the heating process only involves 
the electromagnetic noise of the circuit at frequency Ej, 
while the low frequency noise, which is the dominant 
heating factor in the previous scheme, doesn't affect the 
system. 



DC COOLING 

Effectively, AC driving can be achieved applying a DC 
transport voltage to an auxiliary Josephson junction. We 
modify the system as shown in Fig. [5] so that it becomes 
effectively an SET transistor. Dissipative Cooper pair 
and quasi-particle transport is similar systems was con- 
sidered in Refs. 0,^3- The Hamiltonian in the charge 
basis reads 



H = 



(Q - C R V - C g Vg - C x {x)V x f 
2C S 

cos 6 - Ej, K cos(6 + 2eVt/h) + H x , (37) 



where 6 is the phase on the island (and simultaneously 
on the left junction as the left lead is grounded), Cs = 




FIG. 9: The system with DC driving. 



Cl + Cr, + C g + C x (x). Due to the transport voltage V, 
the Hamiltonian is time dependent and, thus, represents 
a driven system. We repeat the steps described above 
and arrive (in the charge basis) at 



H 



E 



J.L 



X 



a z + hujoa'a + A (cr + a) a z 
(a+ e^-W + hx .) f ( 38 ) 



2 2 



where X = (C R 5V + C B SV e + C X 6V X ), 6<f> = 

(2e/K) J SV dt, and the Josephson frequency uij = 
2eV/H. We see that the right junction's Josephson energy 
assumes the role of pumping amplitude f2, while wj is the 
pumping frequency. Unfortunately V is noisy and, thus, 
the pumping source has a substantial line-width. The 
transport voltage should not be sensitive to the back- 
ground charges, therefore V is assumed to have Ohmic 
noise spectrum (no 1/f component). 

After a 7r/2 rotation in the x — z plane, a x <-> o~ z , we 
obtain 



H 



<J Z - — a x + HuJoa^a + A (a 1 + a) a x 



{{O-Z + i<Jy) e *"*+** +h.C.) 



(39) 



While the Hamiltonians l|12|) and l|39|l look similar, 
there are two important differences. One, already dis- 
cussed, is the fact that the pumping frequency wj in 
(|39|l is fundamentally noisy, while uid in l|12|) can be 
made coherent. The second (very important) difference 
is that in i|12|) the pumping is applied to a z only, while 
in 139fl it couples to a z and a y . Both these facts hin- 
der the cooling. Indeed, the coupling to a v gives a di- 
rect matrix element £j,r/4 between the states |t,n., m) 
and ||, n, rn + 1). This interaction repels the levels and 
we must choose Ej.r -C 4cjo so that the resonant de- 
tuning as in Fig. |21 is possible. In addition, the noise 
of the transport voltage translates into the line width 
for the transition |f,n, to) — > |J_, n — 1, to + 1) equal to 
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= 27i atrkftT '/%, where at r = R/Rq. The fluctuations 
of the transport voltage are not screened by the ratio of 
capacitances as it happens for the gate charge. There- 
fore a tr « 10~ 2 . Because of these additional constraints 
the applicability of the DC cooling scheme is limited to 
higher frequency/quality factor resonators. For an esti- 
mate, consider an oscillator with ujq = 2tt x 1 GHz « 
5 jtteV « 50 mK at temperature T = 50 mK. We then 
obtain T v « 0.3 ^eV, which significantly exceeds T r . 
Hence, we have to substitute T r by T v in all formu- 
las. For the Josephson coupling in the right junction 
we take Ej b, = 2 fieV. Then, instead of Eq. I|15fl. we find 
A w £7 JiR A/(2Sj,l) « 2 • 10~ 3 ^eV (we assume £7j, L w 
50 /LteV). The cooling rate can again be represented as 
where ^_ w 2A 2 /T V « 2 • 10" 5 ^eV. Thus, cooling 
becomes possible only if Q > a;oA4- ~ 2.5 • 10 5 . 



||,n + l> 
r-n> ^ 




li,n-l> 

FIG. 10: Level structure of the voltage driving scheme for ac 
cooling. The solid lines indicate the cooling process; and the 
dotted lines indicate the heating process. 



DISCUSSION AND COMPARISON WITH 
QUANTUM OPTICS. 

For comparison with quantum optics cooling schemes, 
we present here an analysis of some of the processes de- 
scribed above using the quantum optics language. In 
quantum optics literature one usually employs the trans- 
formation to the rotating frame and/or other canon- 
ical transformations together with the Rotating Wave 
Approximation (RWA) in order to single out the near 
resonant terms responsible for the studied transitions. 
We start with Hamiltonian (|12|l and transform it into 
the interaction representation with respect to Ho(t) = 
— [Ej + D, + O cos 2^d t] a z /2 + uioa^a. This amounts to 

a+ -> CT+e -*(*J*+n*+a£ 8inaudt > and a ^ ae — o* . 

(40) 

As il <C 2cjd ~ Ej we can expand the factor 
exp(— i(fi/2w<i) sm2w&£). The near resonant (RWA) 
parts of the coupling term H\ = \(a' + a)a x read, then, 
after the transformation as 

H RWA = {a(T _ + a)(J _ e 2^ t + h c ) ^ (41) 

where the resonance condition 2ojd = Ej + Q — ljq was 
assumed. The first term of (|41f) clearly corresponds to 
the resonant transition provided by the matrix element 
Ay/n (see Eq. l(T5)ll. The second term corresponds to the 
off- resonant transition. In analogy with atom optics Q, 
we can introduce Lamb-Dicke parameter rj = X/2ujd ~ 
X/Ej. Note that, unlike the case of trapped atoms/ions 
we do not have near resonant terms flipping just the spin, 
~ U,a± exp(±ZLJo)- Such an interaction would repel the 
levels, and, for a resonant detuning, one would need f2 <SC 
luq. In our case, a much stronger driving is allowed. We 
do, however, encounter this limitation in the DC cooling 
scheme due to the presence of the near-resonant term in 
the Hamiltonian Eq. i|39|l which is proportional to a y . 



Also, in the dissipative term Hx = — (X/2)a x , we find 
a slow contribution: 

= ( ff - e " 0t + h - C ') ' (42) 

which is responsible for the heating process l|18(l dom- 
inated by the low frequency noise. In quantum optics 
this contribution is typically neglected as there are no 
strong low frequency sources. 

In the case of ac voltage cooling, in the rotating wave 
approximation the Hamiltonian is 

ff = —a z + Haj aa 

+ — (a) a + + aer_ + a) O- + ao + ) (43) 

__ ( a+ e- iuJdt + as^) 

where the interaction includes the resonant coupling 
a'<J+ + acr_ between the states ||, n) and ||, n — 1), and 
the off resonant coupling a^er_ + aa + between the states 
||,n) and |J.,n+l) with an energy difference of 2hajQ. 
The time dependence in the last term shows that only 
high frequency fluctuations on the order of Ej induces 
relaxation. 

The voltage driving scheme for ac cooling presents a 
direct analogue to the laser cooling of quantum optics 
systems as is obvious from Fig. 1101 however, it is a sim- 
pler scheme as the coupling is a direct bilinear coupling 
instead of the polaron coupling in quantum optics. 

"DOPPLER" COOLING 

For smaller oscillator frequencies loq and/or larger 
qubit relaxation rates r r , it may happen that the vi- 
brational sidebands are no longer resolved, i.e., T r > luq. 
Then, cooling to the ground state is impossible but the 
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less effective "Doppler" cooling is still feasible. In the 
standard Doppler cooling, cold atoms moving in a laser 
field experience different light pressure depending on 
their velocity. For instance, if the laser is red-detuned 
relative to the atomic transition, then atom moving to- 
wards the light source will be absorbing photons more 
than an atom moving away from the source, which will 
result in a velocity-dependent force on the atom. By 
properly arranging multiple lasers it is then possible to 
achieve significant reduction of the atom's effective tem- 
perature. The same reasoning applies to an atom in 
a trap, when the trap frequency is much smaller than 
the atomic transition line width. Equivalently, in this 
case, the Doppler cooling can be reinterpreted in terms 
of the cooling and heating transition rates for the oscil- 
lating atom. For red-detuned laser, photon absorption 
processes with simultaneous phonon emission dominate 
the ones where an additional phonon is created. One can 
show that, similar to the free atom case, this leads to 
cooling down to temperatures proportional to the atomic 
transition line width. In this section we will demonstrate 
that an analog of Doppler cooling can be performed in 
both AC and DC cooling setups described in the previous 
sections in the regime of non-resolved sidebands. 

AC cooling scheme. We start with the AC-driving 
scheme described by Hamiltonian (jl2(l . In the Doppler 
case it is enough to study, first, the spin's dynamics with- 
out the oscillator. We perform the following transforma- 
tion 



<t + — > <T-|_e 



-z(2ajd t+ ^j— sin 2u) d t) 



(44) 



and keep only the RWA terms in the coherent part of the 
Hamiltonian but all the terms in the part describing the 
interaction with the bath. The transformed Hamiltonian 
reads 



RWA 



8u> 



X 



-7rcr z - — (a+e 



— liu^t 



. * n x 

+ h.c.) + - — — o~ x , 



4w d 2 



(45) 



where Slo is the detuning 5u> = 2wd — Ej — Q. The last 
term in the Hamiltonian, as before (see Eqs. i|19fl an d 
(I18J0 . generates transitions between the ground and ex- 
cited states of the qubit. The corresponding rate, how- 
ever, is typically much smaller than the relaxation rate 
generated by the second term of the Hamiltonian, and 
hence can be neglected. This is in contrast to the AC 
ground states cooling scheme discussed above, where the 
qubit heating rate determined the lowest possible occu- 
pation number of the resonator. The slow (RWA) part 
of the interaction Hamiltonian reads 



H 



ilow 



W t 

-; — A (a 1 -fa) o x 



(46) 



With Hamiltonian l|45|l we write down the Bloch-Redfield 
equations [Hi , ll^j for the spin's density matrix p = Rp, 



where R is the Redfield tensor. In this equation, the 
density matrix is treated as a four-vector. We choose the 
representation p = (l/2)soo + 7c 2 + ao+ + a*cr_, which 
gives p = (s, 7, a, a*) and Trp = s. For proper density 
matrices s — 1. In this representation 



R = 



( 

2 1 r 








-iSuj - 



2 







iSuj —^f) 



(47) 



\ o o o 

The rate T r originates from the second term of (|45|) . 



r r = 



U! — 2u)d I 



(48) 



(cf. Eq. 0). For simplicity we have assumed that at high 
frequency 2wd the temperature is effectively zero, i.e., 

Next we employ the "quantum regression theorem" 
(approximation) to obtain (for t > 0) the correlator 



C x {t) = (a x (t)a x (0)) = Ti(a x e 



cr x Pa 



(49) 



where p^ is the stationary density matrix, Rpoo = 0. 
For t < we can use C x (—t) — C*(t). In the four- vector 
representation the operator a x multiplying from the left 
is given by 



/ 


1/2 -1 

yi/2 i 



1 1 \ 

-1/2 1/2 



/ 



(50) 



Finally, after the Fourier transform we obtain 
C x {w) = -2ReTr (a x [R + iu]^ 1 a x p^) , 



(51) 



which is easy to calculate (using Mathematica). Expand- 
ing near u) = we obtain C x (u>) = S x + tj x u>, where 



S x 



4r r 



and 



48iu 2 + r? 



32SljT t 



(52) 



(53) 



From the spin correlation functions, we can now de- 
termine the transition rates for the resonator, A± = 
(HX/4uj) 2 C x (^fuJo)- Note, that the same expressions for 
A± can be obtained from Eq. i|ltjfl generalized to arbi- 
trary detuning. Then, the secondary "heat bath" intro- 
duced by cooling is characterized by the effective tem- 
perature 



S x 
2r] x 



4Suj 2 



16<5w 



(54) 
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Optimizing with respect to the detuning Slo, we find 
that the minimum (positive) temperature T m - m m T r /4 
is reached for the optimal red detuning SLo opt fts —T T /2. 
For the effective quality factor due to the spin we obtain 



1 



A, 



lo 



fl 2 X 2 V: 



(55) 



The final resonator occupancy can now be obtained from 
Eq. (|33[1 . Clearly, Q* must be smaller than the oscilla- 
tor's own quality factor in order for cooling to be effec- 
tive. For the same parameters as used in the AC cooling 
section, we find that at the optimal detuning Q* ~ 1, 
which corresponds to a nearly overdamped regime, simi- 
lar to optical molasses in atom optics. Hence, a resonator 
with a frequency ujq < T r /2 — 300 kHz can be cooled 
down to temperature T min w T r /4 ~ 10 [iK. Note, that 
the effect of the dominant qubit heating mechanism due 
to the 1/f noise, which we neglected here, if necessary, 
can be managed by reducing the driving strength (at 
the expense of reduced cooling power, i.e. larger Q*). 

Similar analysis applies to the voltage driving scheme 
for AC cooling in the Doppler regime. For X^/n/2 <C 
IY/2, from Eq. flfy 



X 2 T r 



\ 2 T r 



T2 + 4 (Su - ujoY 



(56) 

where Su = Ud — Ej is the detuning of the driving fre- 
quency from the qubit 's ground-to-excited-state transi- 
tion. Cooling is achieved when Su < 0. Neglecting the 
effect of the finite intrinsic Q-factor, the final phonon 
number from Eq. (|32Jl is 



r 2 . + 4 {Su + u f 
{U)f ~ 16\Su\u 

The scheme provides optimal cooling at Su. 



(57) 

opt — r r /2 

with the temperature T m i n = r r /4, same as the flux driv- 
ing case. 

DC cooling scheme. Analysis of the DC cooling scheme 
proceeds analogously. We start with Hamiltonian l|39|) 
and perform transformation 



(58) 



where, for now, Q = -Ej,r, and we remind that uj = 
(2e/h)V. For simplicity we consider V noiseless now but 
later will introduce the low-frequency (classical) noise of 
V. Then, we obtain 



RWA 



Su 

x t - 



v* + h.c) + SL£<r mt (59) 



where Slo = uj — Ej^. The main differences with <|45[) 
are: i) the second term of l|59|l is absent in 145|l (recall 



the extra matrix element in the DC scheme); ii) the de- 
tuning Slo is noisy due to the noise of V. This will give 
an extra "pure" dephasing rate T v — 2ira tr k-BT '/h, where 
a tI = R/Rq. As in the AC Doppler case, neglecting the 
excitation rate coming from the last term in the Hamil- 
tonian, for the Rcdfield tensor we obtain 



R = 



( 

9 1 r 







ill 

4 





in 

4 



_ m 

2 

2 



-iSu - 




(60) 



iSu — f / 



where T = + T T /2. The general expressions for S x 
and rj x are quite complicated. An analysis shows that the 
simplest cooling regime is achieved when O < T r k T^. 
That is, it does not make sense to keep small T r when the 
line width is anyway given by large T v . The relaxation 
rate T r can easily be increased by choosing bigger gate 
capacitances C x and/or C g . Assuming T r 



T v and O < 



T v we obtain 



and 



12 r„ 



4Su 2 



9T2 



Vx 



96Su T u 



(4Su 2 



9T 2 ) 2 



(61) 



(62) 



Optimizing with respect to Slo, we find that the minimum 
temperature T m ; n « (3/4) T v is achieved for the optimal 
detuning Slo w — (3/2)T v . The effective quality factor 
due to the cooling environment, similarly to i|55|) . is 



n 2 \ 2 Vx 
2u 2 



(63) 



The regime r r <C T v is more subtle and requires further 
analysis. 

Example: Choosing C x » 200 aF instead of previously 
assumed 20 aF, we obtain r r w 0.3 /xeV « r», We choose 
fi = £j.r ~ 0.1 /xeV (note that this is quite a small value 
for usual Josephson junctions). Then we obtain T m { n w 
0.25 /xeV w 2.5 mK for Slo w 0.5 ^eV « 2vr x 100 MHz. 
Near optimal detuning point, we obtain rj x « 6 (/ieV)~ 2 . 
For the coupling constant A we can take A w 1 fj,eV in- 
stead of 0.1 /ieV as we have allowed ten times bigger ca- 
pacitance C x . Then we obtain Q* ss 0.8 • 10 5 . Thus, for 
cooling to be effective, resonator Q-factor should exceed 
10 5 and the final resonator temperature is determined 
according to Eq. (J33J) with T* = 2.5 [iK. 



CONCLUSIONS 

We considered several approaches to active cooling of 
mechanical resonators using a coupling to a supercon- 
ducting Josephson qubit. In the resolved vibrational 
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sideband regime, when the qubit level width is smaller 
than the resonator frequency, we proposed two schemes 
for ground-state cooling of the resonator. In the first 
scheme, the AC driving required for cooling is provided 
by an external microwave source. We find that for a 100 
MHz oscillator coupled to a practically realizable Joseph- 
son qubit, at the external temperatures below 1 K, it is 
possible to reduce the thermal occupancy of the oscillator 
mode by three order of magnitude. In the second scheme, 
the AC driving is generated by the AC Josephson oscilla- 
tions on an auxiliary junction of the qubit. This scheme 
is attractive since there is no need for an external AC 
driving source; however, in the present realization, we 
find that it is not as effective as the one with an explic- 
itly applied AC driving. We also demonstrate that even 
in the regime when the vibrational sidebands are not re- 
solved, it is possible to perform an analogue of Dopplcr 
cooling with the final resonator temperature limited by 
the qubit line width. 
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